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Abstract  (20.) 

Classified  according  to  their  Target  Seeking  or  Orbital  characteristics. 
Flight  paths  in  the  Target  Seeking  domain  are  shown  to  be  convergent 
when  the  system  has  a  wind  penetration  capability.  The  family  of 
Target  Orbiting  solutions  tire  shown  to  produce  trajectories  which  are 
captured  about  the  target  in  Elliptic  spirals.  Launch  criteria  are 
established  from  tha  time  solution  which  takes  the  form  of  an  ellipse 
for  all  the  categories  of  Azimuth  Homing.  The  definition  of  a  Release 
Path  as  the  locus  of  points  from  which  the  flight  time  necessary  vo 
reach  the  target  is  constant,  follows  from  this  result. 
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FOREWARD 


This  report  presents  an  analytic  treatment  concerned  with  a  particular  guidance 
technique  applicable  to  managing  the  flight  path  of  gliding  vehicles  or  gliding  decelerators. 

This  study  was  conducted  under  Department  of  the  Army  Project  No. 
1F2662203AA33.  Guidance  and  Control  of  Gliding  Decelerators. 
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NOMENCLATURE 
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L/D 
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V 
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r 


X 

Subscript: 


Lift  to  Drag  Ratio 
Time 

Horizontal  Component  of  Total  Airspeed  Vector 

|  U  |  ,  Horizontal  Airspeed 

System  Total  Airspeed  Vector 

Radial  Offset  Angle 

Wind  or  Field  Velocity 

|  W  |  ,  Wind  or  Field  Speed 

Horizontal  Space  Coordinate  Fixed  to  Earth 

Horizontal  Space  Coordinate  Perpendicular  to  x  and  fixed  to  earth 

Vertical  Space  Coordinate  Perpendicular  to  the  x-y  Plane 

Magnituda  of  the  radius  Vector  in  Polar  Coordinates 

Angular  Position  in  Pola:  Coordinates 

u/w  Wind  Penetration  Parameter. 

Initial 
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INTRODUCTION 


Automatic  control  of  a  vehicle  in  gliding  flight  poses  some  interesting  problems 
particularly  when  aerodynamic  performance  is  limited.  A  remotely  guided  recovery  system, 
utilizing  some  form  of  gliding  deceleration1,  is  an  example  of  such  a  case.  Typically, 
these  systems  operate  at  relatively  low  airspeed  with  virtually  no  capacity  for  L/D 
modulation.  Discrete  or  continuous  regulation  of  flight  direction  then  becomes  the 
principal  means  of  trajectory  control.  Consequently,  investigations  concerned  with  this 
specific  aspect  of  flight  path  management,  take  on  practical  significance  with  respect  to 
evaluating  the  capabilities  of  various  steering  or  homing  techniques  applicable  to  the 
guidance  of  an  unpowered  gliding  vehicle. 

Analysis  dealing  with  controlled  gliding  flight,  can  cover  a  broaJ  spectrum  ranging 
from  the  application  of  Optimal  Control  Theory1  to  investigations  concerned  simply  with 
the  geometry  of  the  motion.3  Contrasting  the  numerical  trajectory  determinations  which 
follow  from  treatments  such  as  those  developed  in  References  2  and  3,  the  results  of 
this  study  are  analytic  in  scope.  A  class  of  closed  form  solutions  are  derived  from  kinematic 
considerations  of  a  particle  maneuvering  with  constant  speed  through  a  uniform  velocity 
field.  This  representation  serves  to  approximate  the  behavior  of  a  gliding  system  executing 
moderate  turns  in  a  constant  wind  environment.  Assuming  a  specific  guidance  law,  or 
what  amounts  to  a  directional  constraint  on  the  airspeed  vector  eliminates  dynamic 
considerations.  Consequently,  velocity  relationships  can  be  derived  directly  and  solutions 
for  position  and  time  coordinates  obtained  by  integration.  These  results  serve  to  quantify 
the  performance  of  the  assumed  guidance  law  while  providing  valuable  insight  into 
generalized  capabilities. 


ANALYSIS 

The  guidance  law  or  controller  used  in  this  formulation  is  thought  of  as  one  which 
causes  the  system  to  maintain  a  fixed  angular  orientation  between  the  horizontal  projection 
of  its  airspeed  vector  and  a  radial  line  connecting  it  to  the  intended  target.  Figure  (1) 
depicts  the  planar  geometry  of  this  motion,  introducing  the  idea  of  what  shall  be  now 
termed  "Azimuth  Homing".  The  system  is  idealized  as  a  point  mass  and  under  the 
combined  assumptions  of  constant  L/D  and  moderate  turns,  only  the  horizontal  coordinates 
remain  coupled.  The  vertical  mode,  not  shown  in  Figure  (1),  Is  linear  with  time  t,  as 


Indicated  by  equation  (3).  The  vector  quantity  U  defines  the  horizontal  component  of 
the  total  airspeed  vector  V.  U  is  constant  In  mapitude,  however,  Its  direction  may  be 
varied  according  to  control  inputs.  Azimuth  Homing,  as  defined,  constrains  U  to 
continuously  point  along  the  current  radial  at  some  fixed  angular  offset  a.  The  wind, 
or  field  velocity,  W,  is  assumed  to  be  steady,  to  lie  entirely  in  the  x  -  y  plane,  and 
to  point  along  the  y  axis  in  a  negative  sente.  Motion  through  the  wind  field  serves  to 
upset  kinematic  equilibrium  and  as  a  consequence,  the  system  is  always  in  a  state  of 
maneuver. 

It  is  now  possible  to  derive  a  vector  equation  relating  the  absolute  velocity  of  the 
system  relative  to  an  earth  fixed  reference,  to  the  sum  of  V  and  W.  Expressed  in  cylindrical 
coordinates  the  scalar  equations  obtained  from  this  vector  equality  are: 


dr/dt  =  -w(Xcosa+sin0), 

(1) 

(r)d0/dt  *  -w(Xsina+cos0), 

(2) 

dz/dt  ■  (— u)/(L/D), 

(3) 

where  r,  0,  and  z  are  the  conventional  cylindrical  coordinate  designations. 

TRAJECTORY  DETERMINATION 

General  Formulation:  Separating  the  time  dependency  from  equations  (1)  and  (2)  yields 
the  expression; 

dr /r  ■  (Xcosa)d0/(Xsinof+cosfl)  -d(cosfl)/{Asinor+cos0).  (4) 

Motivated  by  physical  considerations,  solutions  to  this  differential  form  may  be  divided 
into  two  categories,  classified  according  to  their  "Target  Seeking",  or  "Target  Orbiting" 
properties.  These  classifications  occur  naturally  reflecting  the  system's  capacity  to  sustain 
angular  motion  with  respect  to  the  target.  This  ability  Is  evident  from  equation  (2)  through 
investigation  of  the  parameter  Xsina.  When  |  Asina  |  <1,  conveigence  to  a  particular 
ray  results.  Given  sufficient  effective  wind  penetration,  the  system  will  seek  this  stable 
angular  alignment  as  r  approaches  zero.  Hence,  the  terminology,  'Target  Seeking". 
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Alternatively,  cyclic  or  orbital  paths  can  be  described  when  |  Xsina  |  >1.  This  velocity 
condition  allows  for  continuous  angular  motion  about  tha  homing  point  from  all  locations 
in  the  wind  field. 

Prior  to  integrating  (4),  a  simplified  and  physically  revealing  version  can  be  obtained 
through  the  variable  transformation; 

cost?  =  {cos 0  -  e)/(1-ecos0),  (5) 

where  0  is  the  so  called  "Eccentric  Anomaly",  a  nomenclature  originating  from  the 
development  of  this  equation  in  Orbital  Mechanics.  The  magnitude  of  e  defined  by  (5) 
is  constant  for  values  between  zero  and  one,  and  for  the  purposes  of  Azimuth  Homing 
assumes  the  role  of  Xsina  or  (1 /Xsina)  depending  upon  the  application. 

Target  Seeking  Trajectories:  e  =  Xsina  for  |  Xsina  |  <1.  Transforming  (1),  (2)  and  (4) 


to  r,  j3  space  yields; 

dr/dt  =  -w(1-e2)yM[(X2-e2)/(1-e2)]1/i  +sin/3/(1  — ecos^)  ] ,  (6) 

rd/3/dt  =  -w{1-e2)1/4cos0,  (7) 

dr/r  =  [(X2-e2)/(1-e2)] 1/4  (d0/cos0)^(cos0)/cos0(1-ecos0).  (8) 

Equation  (8)  is  integrated  directly  to  give; 

r  =  Ksec0{1— ecos0)[se<$+tan0]  ^2— f2  ^ 

where  K  is  an  integration  constant. 


The  effect  of  the  0  transformation  can  now  be  appreciated  by  allowing  e,  in  equations 
(5)  through  (9)  to  approach  zero  white  X  is  constrained  to  remain  finite.  In  terms  of 
physical  coordinates  this  defines  the  special  case  of  "Radial  Homing".4  A  comparison 
betwean  equations  (6)  and  (9)  for  arbitrary  e,  with  those  produced  when  e  =  0  reveals 
their  similarity,  particularly  regard’nq  equations  (7)  and  (9).  It  can  be  concluded,  based 
on  this  similarity,  that  Azimuth  Homing,  constrained  such  that  |  Xsina  I  <1,  ir  a  general 
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f^rm  of  Radial  Homing  when  viewed  in  r  -  0  space.  Therefore,  the  terminal  state 
capabilities  of  Radial  Homing  derived  from  the  physical  plane  establishes  these 
characteristics  for  an  entire  Azimuth  Homing  family.  The  significant  features  of  thsse 
trajectories  are  summarized  as  follows. 

Launched  at  sufficient  altitude,  an  Azimuth  Homing  system  will  reach  the  target 
or  homing  point  prior  to  impact  provided  it  has  the  ability  to  penetrate  the  wind,  i.e. 
when  the  systems  "Generalized  Wind  Penetration  Parameter"  [Xa  — ea  /I — ea  ]  ^>1.  At 
target  arrival  the  absolute  or  resultant  velocity  vector  will  be  aligned  in  a  negative  sense 
alof  *  he  p  m  -*I7  my.  The  arbitrary  ciignment  of  the  wind  velocity  with  the  y  axis, 
which  assisted  in  the  formulation  of  the  basic  theory,  has  no  bearing  on  these  results. 
Any  other  axis  selection  merely  constitutes  a  rotation  of  the  resulting  trajectory  relative 
to  these  priman'  coordinates  producing  no  effect  on  the  end  state.  The  terminal 
characteristics,  (r*0,  /3=-ir/2),  are  also  essentially  independent  of  initial  homing  conditions 
determined  by  the  algebraic  sign  of  e. 

Orbital  Trajectories:  e  *  (1/Xsina)  for  |  (lAsina)  |  <1.  Transforming  (1),  (2)  and  (4) 


to  r  -  0  space  yield:; 

dr/dt  -  (— w/e)  [  1  -ea  ] % [ b/<  1  - c3  )%  +«in/3/(  1  -ecosfl)  ] ,  (10) 

rd/3/dt  -  (-w/e)[1-c,]%.  (11) 

and,  dr/r  -  (b/[1-ca]%)d^-d(coos0)/(1-ecos0),  (12) 

where  b  ■  etna. 

Equation  (12)  is  integrated  directly  to  give; 

,  -  "31 

where  G  is  an  integration  constant 
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The  net  behavior  of  (13)  can  be  determined  by  considering  separately,  the  properties 
of  the  trigonometric  and  exponential  facton  which  combine  to  form  r.  The  quantity 
(l-ecos0)  is  identified  as  a  polar  representation  describing  one  nf  the  Limacon's  of  Pascal, 
thereby  defining  a  closed  curve  in  the  r  —  0  plane.  The  equivalent  figure  in  physical 
coordinates  is  an  ellipse  and  comprises  the  entire  trajectory  when  a  is  allowed  to  assume 
either  of  its  extreme  values  (i.e.  when  a  *  ±ir/2).  Overall,  angular  motion  relative  to 
the  target  can  be  positive  or  negative  depending  upon  the  arbitrary  s‘gn  of  the  radial 
offset  angle  a.  Regardless  of  the  sense  of  change  in  0  or  0,  the  argument  of  e  remains 
negative.  Consequently,  the  exponential  terms  always  acts  as  a  damping  factor,  continually 
suppressing  radial  excusions  in  successive  cycles  about  the  focal  point.  Given  the  nature 
of  the  functions  defining  (13),  the  trajectories  so  generated  may  be  properly  termed 
"Elliptic  Spirals". 

LAUNCH  CRITERIA 

It  is  now  desirable  to  obtain  solutions  in  terms  of  the  time  or  altitude  coordinate. 
Such  relationships  will  be  necessary  in  order  to  identify  spaciai  locations  compatible  with 
the  trajectories  previously  obtained.  Returning  to  equations  (2),  ucd  (3); 

dt  =  (-L/D)/u)dz  *  [(-1/w)rdfl]/[Xsincr+cose) ,  (14) 

where  r  is  now  a  known  function  of  0,  or  0. 

Solutions  to  (14),  which  exist  at  least  in  principle,  can  now  be  attempted  for  both 
the  Target  Seeking  and  Target  Orbiting  cases. 

Time  Integral: 

Two  forms  of  valuation  (14)  are  generated  by  the  respective  utilization  of  equations 
(7;  and  (9)  or  equations  (11)  and  (13).  The  Target  Seeking  form  of  (14)  is; 

(L/D)  (zj/u)  *  (-K/w[1-eJ  j’^)  [  /  sec1^(1-ecos^)(secj3+tan9)^J_e2^^1_eJ^  ^d0j. 

(15) 

The  Target  Orbiting  Version  of  (14)  becomes; 

(L/D)  (zj/u)  -  -Ge/w[1-eJ]%  /  (1-ecos0)  e{b/^~e2^d0,  (16) 
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where  in  both  ca9es  integration  on  altitude  passes  from  the  initial  point  z;  to  the  ground 
plane  (z= 0).  Radial  and  angular  coordinates  proceed  from  the  initial  state  r;,  6\  to  final 
positions  r,  0  which  now  specify  the  systems  location  at  impact.  Both  (15)  and  (16) 
can  be  integrated  by  parts  producing  the  following  solution  which  holds  for  either 
circumstance. 

(L/D)(zjCoso/u)  =  [A/(A2— 1)1  [X— sin(0j— a)]  [1— (r/rj)(A— sin(0— a))/(A— sin(0j— a))J ,  (17) 

where  the  position  coordinate  r  is  given  by  (9)  for  trajectories  in  the  target  seeking  domain 
or  by  (13)  for  the  target  orbiting  case. 

Release  Path  Definition: 

Given  a  specified  launch  or  initial  altitude  z;,  system  parameters  in  terms  of  L/D, 
a,  and  X  which  reflects  wind  properties  as  well,  the  following  physical  interpretation  can 
be  attached  to  equation  (17).  Allowing  for  sufficient  performance  in  terms  of  either 
wind  penetration  ability,  or  the  potential  to  execute  multiple  orbits  about  the  homing 
point,  r  can  be  made  to  approach  zero  in  (17).  Then,  relative  to  a  known  nominal  wind 
direction,  the  curve  so  described  by  these  equations  defines  the  locus  of  points  (z;,  r j, 
0j)  from  which  an  Azimuth  Homing  trajectory  will  terminate  precisely  on  target.  Equation 
(17),  with  r  *  0  the  final  state,  identifies  this  release  path  or  launch  cun/e  as  an  ellipse. 
It  follows  from  these  developments,  that  trajectories  beginning  inside  the  ellipse  will  arrive 
too  early,  (i.e.  an  excess  altitude  condition),  while  those  initiated  outside  will  fall  short. 

SUMMARY  AND  CONCLUSIONS 

Analysis  of  a  particular  guidance  technique  applicable  to  gliding  vehicle  flight  path 
control  has  been  presented.  The  development  has  proceeded  through  three  phases  as 
follows;  (1)  Formulation  of  kinematic  behavior,  leading  to  the  derivation  of  basic 
relationships,  (2)  Mathematical  transformation  of  variables  with  subsequent  solution  to 
the  governing  differential  equations,  and  finally,  (3)  Classification  and  interpretation  of 
results  according  to  the  dictates  of  key  parameters. 

The  Azimuth  Homing  guidance  technique  has  been  shown  to  possess  desirable 
performance  particularly  regarding  the  ability  to  either  maintain  position  relative  to  a 
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homing  point  in  a  captured  spiral  orbit,  or  to  achieve  absolute  target  convergence  in  the 
horizontal  plane.  Since  guidance  is  derived  exclusively  through  management  of  flight 
direction,  consistent  accuracy  will  not  be  possible  in  general  circumstances.  However, 
a  scheme  employing  this  form  of  guidance  ugmented  by  computations  utilizing  position 
measurements  could  be  used  to  control  the  entire  flight.  The  radial  Offset  Angle  (a), 
being  capable  of  discrete  variation,  can  be  utilized  so  as  to  adjust  the  trajectory  to  account 
for  uncontrollable  changes  in  wind  velocity.  The  Wind  Penetration  Parameter  (X)  may 
also  be  considered  as  a  control  parameter  assuming  some  L/D  modulation  were  possible. 
Theoretically,  and  within  the  special  locations  where  terminal  solutions  to  the  Azimuth 
Homing  problem  exist,  the  entire  landing  state  could  be  controlled. 
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